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Homogeneous incompressible MHD turbulence has been studied in the literature in two 
main regimes: with and without an externally imposed uniform mean magnetic field. In 
the former, explicitly anisotropic, case, Kolmogorov-style phenomenologies H |2|] predict a 
state of detailed scale-by-scale equipartion between hydrodynamic and magnetic fluctua- 
tions, which are Alfven waves propagating along the mean field. With proposed scalings 
for the resulting spectra ranging between k~ 2 and k~ 5 ^ 3 , both the hydrodynamic and the 
magnetic energies are concentrated at the outer (forcing) scale of the turbulence. It has long 
been the commonly held view that the other, isotropic, regime of MHD turbulence would, 
in the fully developed state, be very similar, with the energetic magnetic fluctuations at the 
outer scale providing the effective mean field to support the Alfven waves in the inertial 
range The isotropic regime is closely related to the problem of turbulent MHD dynamo, 
where a weak initial magnetic field is embedded into a turbulent conducting medium and 
its evolution is studied. In view of the arguments above, the dynamo has been expected to 
culminate in an Alfvenic equipartition state. In what follows, we will see this outcome is 
less obvious than it seems. 

Because of the close correlation between velocities and magnetic fields in the Alfven 
regime, the ratio between the fluid viscosity and the magnetic diffusivity (the magnetic 
Prandtl number Pr = vjr\) has not been considered important, as long as both v and rj were 
small enough for the turbulence to develop. However, it is a prominent property of hot low- 
density turbulent astrophysical plasmas (interstellar medium, some accretion discs and jets, 
protogalaxies, galaxy clusters, early Universe etc.) that their Pr tends to be extremely large 
and, therefore, that they possess a broad range of subviscous scales at which magnetic fields 
can exist, while velocities cannot. The magnetic field lines are nearly perfectly frozen into 
such a highly conducting fluid. The fluid motions, even though restricted to scales above the 
viscous scale £ u , can excite magnetic fluctuations at much smaller scales via stretching and 
folding of the field lines. In the kinematic (weak-field) regime, the result is an exponentially 
fast pile-up of magnetic energy at the resistive scale £ v ~ Pr -1 / 2 ^,, with a k 3 ! 2 spectrum 
extending through the subviscous range. The associated time scale is the turnover time of 
the viscous-scale turbulent eddies [|1], |J. 

The salient feature of the kinematically generated small-scale fields is their folding struc- 
ture: the fields are "folded" in such a way that the smallness of their characteristic scale 
is due to rapid transverse spatial oscillation of the field direction, while the field lines re- 
main largely unbent up to the scale of the stretching eddy [H 151 H, 0|. Quantitatively, the 
field structure can be studied in terms of statistics of the field-line curvature K = |b • Vb| 
(where b = B/i?) and of its correlation with the field strength B. In the kinematic case, 
analytic theory is possible subject to certain modelling assumptions [§]. The bulk of the 
PDF of K turns out to be at the velocity scales, with a power tail ~ K~ 13 ^ 7 extending 
through the subviscous range K 3> £~ (this scaling is very well confirmed numerically ||). 
This reflects the predominant straightness of the folded field lines at subviscous scales. They 
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are significantly curved only in the "corners" of the folds — the power tail of the curvature 
PDF describes the intermittent distribution of these corners. Furthermore, the field strength 
is anticorrelated with the field-line curvature, i.e., the field is stronger in the straight parts 
of the folds than in the corners || |!| . 

The most important implication of these results is that the system is in a reduced-tension 
state: the tension force, which controls the nonlinear back reaction, is B • VB ~ KB 2 ~ 
B 2 j£ v during the kinematic stage ||, so the nonlinearity becomes important (B • VB ~ 
u-Vu) when the magnetic energy approaches that of the viscous-scale eddies (in a chaotically 
tangled field, this would happen at much smaller magnetic energies, because we would 
have B • VB ~ B 2 /£ v ). Even after the onset of the nonlinear regime, the field in the corners 
remains weak and cannot be expected to exert a significant amount of back reaction on 
the flow. In fact, the folding structure, once set up, is generally very hard to undo, as 
the detailed "unwinding flows" required for that cannot exist at subviscous scales. These 
arguments are corroborated by the numerical simulations 0, which show the persistence of 
the folding structure in the nonlinear regime. In particular, both the anticorrelation between 
the field strength and curvature and the K~ 13 / 7 subviscous-range scaling of the curvature 
PDF remain unchanged. The rigidity of the folding structure plays a key role in the further 
developments. 

We now propose the following scenario of the nonlinear-dynamo evolution [0. The MHD 
induction equation gives the evolution law for the magnetic energy: 

d t (B 2 ) = 2 1 (t)(B 2 )-2r ] k 2 ms (t)(B 2 ), where ^(t) = (BB : Vu)/ (B 2 ), (1) 

and k Tras (t) is the spectrum-integrated rms wave number of the magnetic field. In Eq. ([I]), 
7(t) can be interpreted as the effective stretching rate at time t. During the kinematic 
stage, j(t) ~ ii£ v /£ u , the turnover rate of the viscous-scale eddies. Once the magnetic 
energy becomes comparable to the energy of these eddies, the Lorentz back reaction must 
act to suppress the stretching motions associated with them. However, the next-larger- 
scale eddies still have energies above that of the magnetic field, though the turnover rate of 
these eddies is smaller. These eddies will continue to amplify the field at this slower rate. 
The scale of these eddies is £ s > £ u , so the folds are elongated accordingly and the tension 
force is B ■ VB ~ B 2 /£ s . The corresponding inertial term is u ■ Vu ~ ujj£ s , so these eddies 
become suppressed when B 2 ~ wf , whereupon it will be the turn of the next-larger eddies to 
provide the dominant stretching action. Thus, at any given time, the "stretching scale" £ s (t) 
is defined by u 2 , t x ~ {B 2 )(t) and j(t) is the turnover rate of the eddies of scale £ 8 (t): 
j(t) ~ u ts ( t )/£ B (t). The first term on the rhs of Eq. (0) is then -f(t){B 2 ){t) ~ u\ {t) /£ s {t) ~ e, 
where e = const is the Kolmogorov energy flux. The physical meaning of this result is as 
follows. The turbulent energy injected at the forcing scale cascades hydrodynamically down 
to the scale £ s (t) where (a finite fraction of) it is diverted into the small-scale magnetic 
fields. We conclude that the magnetic energy should grow linearly with time during this 
stage, (B 2 )(t) ~ et, and j(t) ~ 1/t. We stress that the magnetic field is still organized in 
folds of characteristic length £ s (t) with direction reversals at the resistive scale £ v . Comparing 
the two terms in the rhs of Eq. we can estimate £ v ~ k~^ s {t) ~ ll(t)/ri]~ 1 ~ (r/t) ' , 
so £ rj increases. Indeed, as the stretching slows down, selective decay eliminates the modes 
at the extreme UV end of the magnetic spectrum for which the resistive time is now shorter 
than the stretching time and which, consequently, cannot be sustained anymore. 

Some fluid motions do survive at scales below £ s (t) and down to the viscous cutoff. They 
are Alfven waves that propagate along the folds of direction-reversing magnetic fields and do 
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not amplify the magnetic energy J7|. Their dispersion relation is u = ±feii a/ (-B 2 ), where fcn = 

A A 1/2 

(kk : bb) . A finite fraction of the hydro dynamic energy arriving from the large scales is 
channelled into the turbulence of these waves. Since ^T 1 <C fcy <C u; is larger than the 
resistive-dissipation rate of the small-scale fields: to ^> U£ s /£ s ~ 7 ~ )jfc r 2 ms . Therefore, the 
Alfven waves are mostly dissipated viscously via the MHD turbulent cascade 0, [| , rather 
than resistively. This enables us to consider the evolution of the small-scale fields separately 
from that of the Alfven waves. 

The nonlinear-growth stage continues until the magnetic energy becomes comparable 
to the energy of the outer-scale eddies, (B 2 ) ~ ue , and £ 8 ~ £q ~ Re 3 ^£ u , the forcing 
scale. Thus, energy equipartition between magnetic and velocity fields is achieved. The 
time scale for this to happen is the turnover time of the outer-scale eddies t ~ (m^ /£ ) _1 ~ 
~Re 1 ^ 2 (ue l/ /£ u )~ 1 . At this point, 7 ~ ue /£o- Therefore, the resistive scale is now £ v ~ 
(Re Pr) -1 ^ 2 £ , which is larger than its kinematic value ~ Pr~ 1 ^ 2 £ u by a factor of Re 1//4 . 
As there are no scales in the system larger than £ Q , there can be no further growth of 
the magnetic energy. Stretching and bending the ever more rigid magnetic field becomes 
increasingly harder and, instead of amplifying the magnetic energy, causes the field to spring 
back. Therefore, the rate of the energy transfer into the small-scale magnetic field decreases, 
as this channel of energy dissipation becomes inefficient. Instead, we expect the energy 
injected by the forcing to be increasingly diverted into the Alfvenic turbulence that is left 
throughout the inertial range in the wake of the suppression of the stretching motions. 
During this last stage of the nonlinear dynamo, the magnetic energy stays approximately 
constant (growing only very slightly) while j(t) drops below the turnover rate of the outer- 
scale eddies. The energy balance Eq. ([]]) then implies further movement of £ v ~ k~^ s (t) 
towards larger scales. The mechanism for this decrease is the same as in the nonlinear- 
growth stage: the resistive decay of the high-A; modes outpaces their regeneration by the 
weakened stretching. 

Thus, £ v increases both in the nonlinear-growth stage and during the subsequent slower 
approach to saturation. In [0], we show that the evolution of the magnetic-energy spectrum 
in both regimes is likely to be self-similar with k rras (t) ~ (r/t) -1 / 2 . It is very important 
to understand how far it can proceed. In our arguments so far, we have disregarded the 
Alfvenic component of the turbulence. It is not, however, justified to do so once the small- 
scale magnetic energy reaches the velocity scales. Indeed, the decrease of k rms is basically 
a consequence of the balance between the field-amplification and the resistive-dissipation 
terms in Eq. ([]]). However, once k rms ~ t~ x , the Alfvenic turbulence will start to affect k rms 
in an essential way: since the waves are damped at the viscous scale, fc rms must stabilize at 
fc rms ~ l~ x . The resulting turbulent state features folded magnetic fields reversing directions 
at the viscous scale plus Alfven waves in the inertial range propagating along the folds. There 
are then two possibilities: either (i) this represents the final steady state of the isotropic 
MHD turbulence, or (ii) further evolution will lead to unwinding of the folds and continued 
energy transfer to larger scales, so that the spectrum will eventually peak at the outer scale 
and have an Alfvenic power tail extending through the inertial range. It is in the latter case 
that the turbulence in the inertial range would be of the usual Alfven- wave kind [|3], |2j , where 
the large-scale magnetic fluctuations would provide a mean field along which the inertial- 
range Alfven waves propagate. However, in order to achieve such a state, the folds would 
have to be unwound. In view of the rigidity of the folding structure, it is unclear how this 
can be done. This dichotomy remains unresolved and requires further study. 

In either case, k Tras ~ so the fully developed isotropic MHD turbulence is character- 
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ized by the equalization of the resistive and viscous scales. The time scale at which such an 
equalization is brought about is the resistive time associated with the viscous scale of the tur- 
bulence: t v (£ u ) ~ £l/rj. We immediately notice that, in order for the nonlinear-growth stage 
to run its full course up to the energy equipartion (B 2 ) ~ uj Q , this time scale must be longer 
than the turnover time of the outer-scale eddies: t v (£ u ) ^> t , which requires Pr > Re 1/2 and 
is also equivalent to the condition that £ v <C £ u at the end of the nonlinear-growth stage. 
This constitutes the "true large-Pr regime," which is the one relevant for astrophysical plas- 
mas. In this regime, the magnetic energy saturates at the equipartition level, (B 2 ) ~ u 2 Q . 
From Eq. ([I]), we get an estimate of the amount of turbulent power that is still dissipated 
resisively: ■y{B 2 ) ~ {r]/£ 2 )u 2 o ~ eRe^/Pr <C e, i.e., most of the injected power now goes 
into the viscously-dissipated Alfvenic motions. The other possibility is Pr < Re 1//2 . In this 
case, the self-similar dynamo evolution is curtailed during the nonlinear-growth stage with 
magnetic energy still at a subequipartition value, (B 2 )/uj o ~ t v (£ u )/to ~ Pr/Re 1 / 2 , and 
resistivity continuing to take a significant part in the dissipation of the turbulent energy. 
If no further evolution takes place, this gives an estimate of the saturation energy of the 
magnetic component of the turbulence. In this regime, Pr is not large enough to capture all 
of the physics of the large-Pr dynamo. Most of the extant numerical simulations appear to 
be in this regime (see references in [[J]). 

Thus, we have identified two possibilities for the long-time behaviour of the isotropic 
MHD turbulence: saturation in the usual Alfvenic state || [| and saturation with the 
magnetic energy tied up in the viscous-scale fields. Which one is realized depends on the 
way the small-scale folded fields interact with the inertial-range Alfvenic turbulence. We 
stress that there is no numerical evidence available at present that would confirm that, the 
isotropic MHD turbulence without externally imposed mean field — at any Prandtl number 
- attains the Alfvenic state of scale-by- scale equipartition envisioned in the commonly 
accepted phenomenologies 0, 0. In fact, medium-resolution simulations rather seem to 
support the final states with small-scale energy concentration even for Pr = 1. This does not 
mean that the Alfvenic picture is incorrect per se. However, all existing phenomenologies 
of the Alfvenic turbulence depend on the assumption || that the strongest magnetic fields 
are those at the outer scale. This is authomatically satisfied if a finite mean field is imposed 
externally. However, it remains to be seen if such a distribution of energy is set up self- 
consistently when the turbulence is isotropic. 
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